I. INTRODUCTION 



The propagation of electromagnetic waves through plasmas has been extensively studied from a classical viewpoint. 
This is usually justified in most cases by the physical conditions regarding the plasma temperature and density at 
least when the considered wavelengths are not too short. 

In the classical approach the plasma operates as a sort of active filter both absorbing and distorting electromagnetic 
waves in various ways depending on the plasma parameters, its homogeneity, the presence of an external magnetic 
field and of course the frequency of the wave. To have an idea of this complicated methodology and of the kind of 
results one obtains it is useful to have a look to ref. [[[). It is remarkable that, according to the different conditions, 
various frequency shifts arise. 

Another classical approach, particularly fit for optical, or shorter, wavelengths is that which considers the interaction 
of the electromagnetic wave with a fluctuating medium as a scattering by inhomogcneitics inside it. This analysis 
has been made in refs. ||], ||, Q and leads again to frequency shifts of the incoming radiation depending now also 
on the scattering angle. In particular a blue shift is found for forward scattering, as a manifestation of the Rayleigh 
scattering. 

In principle however the quantum aspects of the interaction of the wave with the plasma should not be overlooked. 
This is likely to be especially true in some situations of astrophysical or even cosmological interest, where the density 
C$ ' of the plasma is rather low and the distances are so high to allow also small effects to pile up. 

In this paper we shall precisely investigate the quantum effects using standard quantum electrodynamics and a 
perturbative treatment, as outlined in section ||. The method will produce also the known classical features of the 
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propagation of electromagnetic waves through a plasma, as we show in section [V . On the quantum side, we take into 



We consider a situation where an electromagnetic wave, plane and monochromatic, coexists with a plasma of 



H . account in section [II the many body nature of a plasma that entrains the appropriate fermionic statistic. 

In the non relativistic limit we find, for a low density locally homogeneous (i.e. homogeneous on the scale of a few 
wavelengths) plasma, a blue shift of the photon frequency (this has the same sign as that of the classical results in 
\q • , H , Q , though now the effect has a completely different origin) . 

Finally in section ^ we discuss the validity conditions of our approximations, in particular those concerning the 
f^*) ■ possibility of overlooking relativistic corrections. 

Just to fix ideas and exemplify our low density plasma, we refer, in various parts of the paper, to numerical values 
of the parameters of the order of those valid for the solar corona, i.e. a number density ~ 10 6 electrons/cm 3 and a 
' temperature ~ 10 6 K. 

df II. BASIC ASSUMPTIONS AND OUTLINE OF THE METHOD 
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22 ■ electrons with a numerical density distribution n(x). The propagation is along the x axis. 

The "unperturbed" state is obtained when the coupling between the wave and the plasma is set to zero; consequently 
^ • the energy distribution of the electrons is that of a fermionic plasma of temperature T and density n restrained into 
' a potential well and the wave has a frequency u>. 

Let us now switch on the coupling and, assuming the interaction energy to be small, we determine the shift in the 
energies as a perturbation of the "initial" (i.e. uncoupled) situation. Supposing that the interaction is set up in a 
finite time lapse and comparing the situations at — oo in time with that at +oo, we mimic the actual process of an 
incoming plane wave of frequency to and an outgoing one of frequency uj + Slu. 

The technique actually used to compute the shift in the frequency of the wave is that of the time independent 
perturbations. The interaction Hamiltonian is: 

Hj = - R - J j • A dr (1) 

where A is the vector potential operator of the wave and j is the current density operator appropriate for this problem; 
the volume integral is limited to the confinement region of the plasma. The complete non relativistic expression for j 
when an electromagnetic interaction is present is g: 

j = -^($ t V$-$V$ t )-— $t$A (2) 
2m mc 

Cast in the second quantization formalism, the interaction Hamiltonian is 

H I = H' I + Hj (3) 
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where we defined 




k,q'-q (6) 



The 6 f 's and the b's are the bosonic creation and annihilation operators associated to the photons; the a's are 
fermionic (electronic) operators. The energy of the photon is hu>, its momentum is Kk; of course lu = ck; the 
polarization of the photon is expressed by the unitary vector u k ; k's and q's are respectively photonic and electronic 
wave vectors; primes denote intermediate state variables; finally, the unessential spin and polarization indices have 
been dropped, though an average over them in the final formulae is performed. 

The last term in (|2|) is usually omitted under the explicit, and some times implicit, assumption it gives contributions 
small with respect to those coming from the first. 

The frequency shift of the external photon can be written in the form 

Su = 8uj' + Suj" 

where 5u>' and 5u" are respectively the contributions due to Hj and H" . In this paper we calculate Slo for a low 
density plasma, with a perturbative treatment up to order a — e 2 /he (this is equivalent to keep the first order term 
in Hj and the second in Hj). 



III. THE QUANTUM FREQUENCY SHIFT 

We first consider the contributions to the shift due to H\. 

Let H em and Hp be the electromagnetic field and the plasma free Hamiltonians; Tl em and Hp are respectively the 
bosonic Fock space associated to the photons and the fermionic Fock space associated to the electrons; <I> k € H em 
and G Tip are the one particle electronic and photonic wave functions. 

We calculate the shift AE Uj q of the value of the energy of the state <E> k <8> of H em ®TLp. Since $ k ® * q is an 
eigenvector of H em + Hp we can use the second order perturbation theory (the first order energy shift is zero) . The 
energy shift of the state i is 

where j labels any (normalized) eigenvector of H ern + Hp and Vj is the probability that the state j is occupied by an 
electron. 

The only non vanishing contributions are those with respect to j states of the form j\ := $o <8> \&q' , j% '■= 
^-^$ k <g> $ k '^ ® ^q' or j 3 := (^-±=(& k , (gi (gi vj/ q/j ($ i s the electromagnetic vacuum). It is easy to show that 

(Hi) hi = flk, 

( H i) h i = (H'l) hi = ^fq',q 

So the energy shift induced by Hj on $ k <8 ^ q is 

AE' = A£« + AE^ (8) 
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Huj + e q — e q > hu> + e q i - e q 
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where, for a Fermi gas, in low density approximation: 
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the electrons ^q. We have 



In order to obtain the energy shift, we have to take the mean value of AE' W with respect to the possible states of 

AE' U = AE^ + AE^ 
AE<?>:=Yv q AEW (12) 



AE^ :=2>,A2i£ 
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We shall consider the contribution to the shift AE^ in the next section. 

The term AE^ gives a divergent contribution |7|. In the following we shall show that, if we subtract the second 
order self energy of the electrons AE° , the contribution to the resulting term is finite and independent from the 
normalization volume. This self energy is 
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Thus, the observable contribution to the shift is: 

AE^ = - AE° 

Using © and (0), 



k' qq' 



e q — € q i — tick' 



(15) 



We now show that y AE^ has a finite value, independent from the normalization volume . 

Since in low density conditions v q v q i = z 2 exp [—j3 (e g / + e q )], we have, using the explicit expression (Q) for the g^, , 
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Summing over q' and using ^ — + * ri J dkj dq, we get 

qk ^ 



1 a 9 e 2 7r 

—AE { 3 ] = z — 5 r 

V w m 2 c2(27r) £ 



dq / dk 



(q ■ u k ) 2 exp (e| q _ k | + e g )] 
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(16) 
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The integral over k can be computed in spherical coordinates (r, 9, <p) choosing q as polar axis; the polarization 
vector (normal to the polar axis) is defined by ip = 0. We have (fi — cos 9): 



(q • u k ) z = q 2 cos 2 ipsva 2 9 = q 2 cos 2 ip (l - fj, 2> 



h 2 

£| q -k| ± e g = w- {o 2 + k2 - 2 1 k V ± I 2 ) 



Calculating the integral ([If]) over ip, then integrating over all the possible directions of q: 

1 ~ m z 2 e 2 h 2 f +oc f + °° , f 1 , ex P r-(2q 2 + k 2 - 2qkfi) Ik 2 ] , s 



with k T := y'S^ . 

In order to calculate the integrals over q and k , we notice that the integrand is not exponentially zero only when 
the positive definite quantity 2q 2 + k 2 — 2qk[i is smaller then a few fc^'s. This condition is satisfied only if k and q 
are small with respect to some fcr's- Then, if the plasma temperature is not higher than ~ 10 6 K , we can neglect 



^ (| - qn) in the denominator of (|T|) with respect to 1 (in fact ^ (f - qp) = O (-^k T ) = O ( J ^ ) , and, if 



T~ lO 6 ^, ~ 1CT 2 ). 

Now the integrals are Gaussian, and they can be easily evaluated . This gives 

I „ z 2 g 2 f2 

Finally in the low density assumption (z <C 1) the density n of the plasma is: 



As a consequence, we have 
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— AE 1 - ' = = a (19) 
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If the plasma is not homogenous, this energy shift depends on the position via the electronic plasma thermodynamic 
parameters. As we noticed in the introduction, the dependence on the position is negligible over distances of the order 
of the wavelength. Thus, the observed energy shift is simply the volume average (see section |v| below) of AE^: 

hSwW / dxA&V (x ) = 



IV. CLASSICAL CONTRIBUTIONS TO THE SHIFT 

We now consider the contribution to the shift due to AE^P . In low density conditions we have, after summing with 
respect to q' and taking the continuous limit, 
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AE 1 = ^ Z-z-zrzs I d**-* 3 " (q ' u k) 2 
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The integral can be calculated in the same approximation that leads to (|T^). We obtain 
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Once again, if the plasma is not homogeneous, 
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where TV is the total number of electrons in V. 

Let us now determine the contribution coming from H l of (^|). 

Using again a perturbative treatment we see that the first order term is no longer zero and its magnitude in terms 
of powers of the coupling parameter a is the same as that of (|20|), it cannot then be a priori overlooked. 
The correction we are now studying is: 
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where H'/ is given by eq. (|). Thus 
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In fact J^ q , u g ' is, by definition, the total number N of particles, consequently y is nothing else than the average 
density (n) of the plasma in the given volume. This gives Su>" — e 2 (n) /(2mck), which is precisely the classical 
correction to the dispersion relation due to the presence of a low density homogeneous electron plasma (whose 



constant density is (n)), when the temperature is zero. In fact, when it is 
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, we have (see M) 
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where d. 



p 



is the plasma frequency. For uj p <C kc and <1 ( non relativistic approximation) this gives 
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This is of course true when T ~ 10 6 K, which implies 
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The last term in (|26|) coincides with 5u>" + 6u>^, where Slo^ gives the first order temperature correction. Like all 
the classical terms, it depends on N/V (V is the normalization volume of the electromagnetic field). Hence, when 
V — > oo , such a term gives a vanishing contribution to the frequency shift. 



V. DISCUSSION 



The dependency of the classical terms on N/V rules them out when the measuring apparatus is outside of the 
plasma and the normalization volume for the electromagnetic wave is infinite. This is not the case for the quantum 
contribution (|20|), which is always different from zero: it will be the only contribution to the frequency shift observable 
in astrophysical or even cosmological conditions. An important remark on (^0|) is that the volume over which one 
integrates cannot be the volume of the whole universe; in fact it extends from the source to the receiver (~ from 
— oo to +oo ) along the line of sight, but transversely it should not be more than the distance over which the plasma 
may practically be thought of as infinite for quantum mechanical calculations. This transversal extension is a sort 
of coherence length for the plasma. It should be less than the wavelength of plasmons, the Debye length and the 
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screening length of the plasma. On the other hand its square should of course be much greater than the Compton 
scattering cross section: our approach has indeed nothing to do with individual scattering phenomena. 
Assuming for simplicity a Gaussian distribution along the line of sight, such as 
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n (x) = riQe ^r- 



© gives: 
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L is the transverse "coherence length": it may in fact be used as a phcnomenological parameter. 

In our model relaxation phenomena in the plasma play no role. This is because we consider a steady state situation 
and, furthermore, our plasma is modelled as a reservoir, that can exchange with the photon any amount of energy 
without changing its thermodynamic state. Thus any dynamic process due to the transit of the photon is neglected. Of 
course, this approximation is good only if the plasma is very extended, its density is very low and the electromagnetic 
field is not too strong, but all of these conditions are indeed satisfied in our case. 

The result we found has been obtained under some assumptions that need a careful consideration. First of all, 
the calculation is based on the perturbation theory. Usually this theory treats single electrons interacting with an 
external field or with a bath of photons. Here we are dealing with one photon interacting with a gas of electrons. The 
fermionic nature of electrons comes into play through the factor v^v^i in (|l5|), or simply, when the approximation 
is appropriate, through the fugacity z of the plasma. When the fugacity is small enough, it is easy to see that the 
relevant terms of the perturbative series are weighed not only by a and its powers, but rather by products of powers 
of a and powers of z: this fact may change the relative importance of the different contributions. 

Actually (^0|) is of order az 2 . If relativistic effects had been taken into account, a second order contribution had 

come also from processes such as virtual pair creation. Such a contribution is proportional to az and can be 

overlooked when 



To sum up we conclude that the blue shift found in this paper applies to waves whose frequency satisfies to the 
condition 



For higher frequencies relativistic contributions must be included. 

To give an idea of the numbers, consider that in the physical conditions of vast portions of the solar corona z can 
be as low as 10 -17 ; consequently the upper frequency is ~ 10 11 Hz. 
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